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The local density of states power spectrum of optimally doped Bi2Sr2CaCu2 08+a; (BSCCO) has 
been interpreted in terms of quasiparticle interference peaks corresponding to an "octet" of scattering 
wave vectors connecting k-points where the density of states is maximal. Until now, theoretical 
treatments have not been able to reproduce the experimentally observed weights and widths of these 
"octet" peaks; in particular, the predominance of the dispersing "qi" peak parallel to the Cu-O bond 
directions has remained a mystery. In addition, such theories predict "background" features which 
are not observed experimentally. Here, we show that most of the discrepancies can be resolved when a 
realistic model for the out-of-plane disorder in BSCCO is used. Weak extended potential scatterers, 
which are assumed to represent cation disorder, suppress large-momentum features and broaden 
the low-energy "q7" -peaks, whereas scattering at order parameter variations, possibly caused by 
a dopant-modulated pair interaction around interstitial oxygens, strongly enhances the dispersing 
"qi"-peaks. 

PACS numbers: 74.25.Bt,74.25.Jb,74.40.+k 



I. INTRODUCTION 

For several years, high-resolution scanning tunnelling 
microscopy (STM) experiments [JH0,SlHllHEIl 
ITol Hit \l2L ITa . have been imaging the local electronic 
structure of cuprate superconductors. These experiments 
have discovered resonant defect states Q, H and re- 
vealed the existence of nano-scale inhomogeneities 0, H 
0,0). From a highly disordered spatial tunneling signal, 
it has become standard to use Fourier transform scanning 
tunnelling spectroscopy (FT-STS) EH techniques to ex- 
tract, e.g., the characteristic wavelengths of Friedcl-type 
oscillations EH and, if present, the background charge or- 
der. In a simple metal, the local density of states around 
an impurity varies as ~ cos2fcpr/r 3 , so the wavelength 
of the local density of states (LDOS) "ripples" caused by 
a single impurity is a measure of the Fermi wave vector of 
the pure system EH- In the quasi-2D d-wavc supercon- 
ductor, the spherical shell of maximum FT-STS intensity 
with radius 2k p is replaced by a discrete set of peaks at 
positions q Q , a = 0...7 (see Fig. which connect an 
octet of points of h igh density of states in momentum 
space H 113 El El El Et The fact that in optimally 
doped BSCCO the measured dispersion of these peaks 
with STM bias |3. ll2l[T^ | agrees semi-quantitatively with 
the predictions based on d-wavc BCS-thcory is consid- 
ered as strong evidence for the existence of well-defined 
quasiparticles in the superconducting state of optimally 
doped cuprate superconductors, and for the applicability 
of d-wave BCS-theory in this regime. Recent attempts to 
relate ARPES spectral functions to STM Fourier trans- 
form patterns have met with reasonable success, support- 
ing this conclusion p i I2I l24j . 

At lower doping, the dispersing peaks are still present, 
but additional nondispersing peaks are observed corre- 
sponding to wave vectors along the Cu-0 bonds with 
wavelength close to four lattice spacings 0,^3- Nondis- 



persing peaks have been reported in overdoped samples 
as well, and have been discussed in terms of stripelikc 
charge ordering d EH • The nature and degree of charge 
ordering and the existence of the nondispcrsive peaks in 
these materials is still controversial, however. In this 
work we neglect charge ordering and discuss exclusively 
the Fourier transform LDOS response of the d-wave su- 
perconducting state to disorder. 

Despite the overall success of interpreting the Fourier 
STM patterns in terms of quasiparticle interference, none 
of the existing models has been able to correctly ac- 
count for the observed weights and widths of the dis- 
persive q-sp ace peaks. The pointlike single-impurity the- 
ories HEJEIEHHEI yield peak structures which are 
sharp, and linear or arclike in structure, in contrast to the 
experimental features, which appear as a series of broad, 
roughly round spots. Depending on the assumed form of 
the impurity potential and its strength, the relative inten- 
sity weights of the different octet peaks change, but there 
has been no choice of potential found which adequately 
accounts for the experimentally observed weights. In par- 
ticular, all such calculations appear to drastically under- 
estimate the weight of the dispersing qi-peaks and in 
contrast to experiment predict many "background" fea- 
tures, i.e., strong intensity in regions of g-space other 
than the octet peaks. 

Many-impurity calculations with weak scatterers dis- 
played some of the same features as the single-impurity 
calculations, and showed how increasing disorder even- 
tually swamps the octet peaks by raising the noise floor, 
but were not notably more successful in reproducing the 
experimental FT-STS patterns po| . It was therefore sug- 
gested by Zhu et al. |2y| that it is essential to consider a 
more realistic disorder model for BSCCO, including the 
effect of roughly 0.2% in-plane native defects, which are 
generally considered as strong pointlike scatterers due to 
the resonances they generate near zero bias E3j as we H 
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FIG. 1: Schematic depiction of the first Brillouin zone in 
cuprate superconductors. Green (blue) regions represent 
Ak > (Ak < 0). Red dots are intersections of constant 
quasiparticle energy contours with the Fermi surface. The 
fraction of octet vectors, qi, q4, qs which should in principle 
be visible in the case of a "pointlike" (4-bond) n-scatterer 
are denoted with solid lines, the remaining octet vectors are 
indicated with dashed lines. 



as the effect of out-of-plane defects, which most likely 
act as weak extended scatterers due to poor screening 
within the Cu02-planes. Using this disorder model Zhu 
et al. were able to reproduce the experimental FT-STS 
patterns at low biases of order 10-15meV reasonably well, 
but at higher energies agreement was rather poor. They 
attributed the low-energy behavior to the in-plane na- 
tive defects and interpreted the role of the extended weak 
out-of-plane scatterers as broadening the octet peaks. A 
broadening of small-momentum features has also been 
reported by Dell'Anna et al. [2]} due to mesoscopic in- 
homogeneities. 

In this paper we show that a realistic treatment of 
the out-of-plane disorder plus the consideration of in- 
plane native defects is necessary to understand the key 
features of the experimental Fourier transform patterns. 
There are at least two major sources of out-of-plane 
disorder in this material: (i) non-stoichiometric oxygen 
dopant atoms and (ii) random substitution of Bi at the 
Sr-sites j2^|. Very recently, direct information on both 
types of impurities has been obtained by STM. Kinoda 
et al. |23| succeeded in imaging the excess Bi 3+ at a 
bias of 1.7 eV, finding an areal concentration of about 
3%, consistent with earlier estimates from bulk measure- 
ments [2S). In addition, McElroy et al. [3(| imaged a set 
of atomic scale defects at -0.96 cV, whose concentration 
scales with doping. These defects, which were ultimately 
identified as the interstitial oxygen dopant atoms, were 
shown to correlate strongly and positively with the local 
gap size in the superconducting state. 

We have recently suggested that the above mentioned 
experimental observations can be understood based on 
the simple assumption that the main effect of the inter- 
stitial dopants (i.e. disorder (i)) is to increase the pair 



interaction locally [31| . This assumption yields gap mod- 
ulations which correlate positively with the dopant po- 
sitions, gives rise to particle-hole symmetric spectra and 
suppresses high-energy coherence peaks as observed in 
experiment. The underlying assumption is that each in- 
terstitial oxygen distorts the "cage" of atoms around it, 
thereby modulating the pair interaction locally. Alterna- 
tively, these order parameter modulations could also arise 
from sources other than microscopic pairing disorder, e.g. 
via coupling of ordinary electrostatic impurity potentials 
to phases which c omp ete with d-wave superconductiv- 
ity, see e.g. Ref. [32]. Nevertheless, the phenomeno- 
logical assumption of microscopic pairing disorder has 
already proven successful with respect to many experi- 
mental observations in BSCCO. Scattering by order pa- 
rameter variations has been reported to improve dramat- 
ically the fit to the inelastic tunneling spectroscopy ex- 
periments [3i|. Fang et al. |34| reported very sharp co- 
herence peaks in small gap regions consistent with an An- 
dreev bound state in a region of suppressed order param- 
eter, as discussed in Ref. jU- Very recently, J.-X. Zhu 
explored the effect of dopant modulated Cu-Cu hopping 
within a Bogoliubiov-de Gennes framework |35j . Finally, 
order parameter modulations with a phase twist have 
been found to yield low-energy impurity bound state en- 
ergies and a spatial distribution of the LDOS similar to 
those observed near Zn- impurities in BSCCO-2212 [3(|. 
Thus it is intriguing to explore the consequences of this 
unusual quasiparticle scattering from the order param- 
eter modulations further, in particular its effect on the 
FT-STS patterns. 

Unlike the oxygen dopants, the Bi 3+ defects, i.e. disor- 
der component (ii), most likely act as conventional weak 
and extended potential scatterers due to the different 
charge of the Sr- and Bi-ions and the smaller deformation 
of the lattice itself. The existence of a large concentration 
of extended weak potential scatterers is also supported by 
an analysis of transport properties in BSCCO |37l |. 

The final aspect of the disorder we discuss is the small 
concentration of in-plane native defects, about 0.2% per 
Cu, imag edby STM It is believed that these are Cu 
vacancies in the CuC"2 plane, although no direct proof 
exists. Individually, these defects provide a resonant sig- 
nal centered within 1 meV of the Fermi level similar but 
not identical to Zn. The resonance modulates the LDOS 
primarily below 15 meV or less, and then furthermore at 
energies not far from the gap edge where the resonant 
spectral weight originates. We model these defects as 
strong pointlike unitary potential scatterers, as in Zhu 
et al. [26J , noting that because the concentration and 
resonance energies of these defects are known from ex- 
periment, they do not introduce new model parameters. 

Thus, we believe that a realistic model of the disor- 
der in BSCCO should contain a net out-of-plane disorder 
consisting of extended weak potential scatterers (approx. 
3%, i.e., the concentration of Bi-dcfccts) and extended 
off-diagonal scatterers (approx. 7.5%, i.e., half the dop- 
ing, because each oxygen dopant most likely contributes 
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two holes) plus a concentration of 0.2% pointlike and 
strong in-plane defects. It is the purpose of this paper to 
investigate the effect of these three disorder components 
on the FT-STS patterns. 



II. MODEL 

As a model for the homogeneous system we use the 
usual BCS Hamiltonian given by 



A k c kT c 



H.c. 



(1) 



k.a 



where £k is the quasiparticle dispersion proposed to fit 
the ARPES data |23: £k = ^Lo tnXn(k), where tos = 
0.1305, -0.5951, 0.1636, -0.0519, -0.1117, 0.0510 (cV), 
and Xo-s(fc) = 1, [cosk x + cosk y )/2, cos k x cos k y , 
(cos2k x + cos2fc y )/2, (cos 2k x cos k y + cos2k y cosk x )/2, 
and cos 2/c.;, cos 2/c,,. For d-wave pairing we have, A k = 
2Ao(cosfca; — cosfcy) with Ao = 12mcV. 

In terms of the Nambu spinor = (cL., c_kj.), the 
corresponding bulk Green's function in Matsubara rep- 
resentation is given by 



iU! n T0 + CkT 3 + A k Tl 



(iU! n 



(2) 



where E£ = £k + A k , and To is the 2x2 identity ma- 
trix, whereas (ti, ra, 13) denote the three Pauli matrices. 
In real-space, the perturbation due to the conventional 
point-like potential scatterers in the diagonal T3 channel 
and modulated pairing amplitude in the off-diagonal T\ 
channel can be written as 

W'(r, r') = $ [v( r - r')S(r - r')r 3 + <SA(r, r')ri] 

(3) 

Note that the conventional impurity scattering occurs 
with Pauli matrix T3 in J3J, whereas the order parameter 
modulation, or Andreev scattering term, occurs with t\. 
We will henceforth adopt the convenient shorthand of re- 
ferring to the two types of processes as T3 or t± disorder, 
respectively. 

We use a Yukawa form to model the extended conven- 
tional potential of an out-of-plane defect: 

y p-d/\ 

V(r) = e _°J A - with d=^T^ z , (4) 

where r is an in-plane vector, r z is the distance of the 
defect from the Cu02-plane and the potential is normal- 
ized to Vq in the Cu02-plane directly below the defect. 
An extended Ti-scatterer is modeled by using the same 
functional form in the t\ channel as 



A A r p~ d / x n 
*A(r, r') = e _ r J x d with d = ^(-(r + r')) 2 + r* . 

(5) 

We note that there is no microscopic justification for this 
form, nor any direct connection with a physical screened 



Coulomb potential; this is merely a convenient way to 
introduce a length scale to the order parameter modula- 
tions induced by the dopant atoms. 

In order to determine the resulting LDOS as a function 
of energy and lattice sites, one needs to obtain the full 
Green's function Q(t,t' ,iu> n ) given by the Dyson equa- 
tion 

S(r,r')=S°(r-r')+ £ g°(r-r")H'(r" , t>")Q{t>" , r'). 

r" .v'" 

(6) 

Thus, by calculating 
G u (r,iuj n ) = y — - -5 exp(«k-r), (7) 



k 



the remaining problem is that of a simple matrix inver- 
sion, and the LDOS p(r,ui) can be extracted from the 
imaginary part of the full Green's function, p(r, u>) = 
-±hn.g n (T,T,u)). 



III. RESULTS 

We now examine the effect of a single scatterer on the 
FT-STS patt erns. The philosophy of calculations of this 
type [rH Il8( is simply that the "Friedel oscillations" in 
the d-wave quasiparticle sea around a single impurity will 
contain many of the wavelengths present in the fully dis- 
ordered system. The FTDOS due to the impurity is then 
given by 



-?qr 



(8) 



rGLxL 



where L x L is a square set of L 2 positions at which 
measurements are made, and q = 27r(m, n)/L are vec- 
tors in the associated reciprocal lattice. The power spec- 
trum |39| is simply the 2D image obtained from the am- 
plitude |p(q, uj)\. In the case that several disorder com- 
ponents are present we calculate the power spectrum as 
|p(q, u)\ = I J^. 7ii/5j(q, u>)\, where pi and rij denote the 
FTDOS and concentration of each disorder component. 

Fig. [5] compares the resulting power spectra for dif- 
ferent impurity models computed for L = 53. In the 
1st column, the power spectrum for a weak pointlike T3- 
impurity, discussed previously in Refs. 0,Hi|, is repro- 
duced. In contrast to the experimental power spectrum, 
the calculation for a pointlike r 3 -impurity predicts strong 
weight at large wavevectors especially near (n, 7r), which 
is observed experimentally only at very low frequencies. 
In the case of an extended T3-scatterer (2nd column in 
Fig. [3), on the other hand, the high- momentum features 
are suppressed due to the rapid fall-off of the impurity 
potential V(q) in momentum space. This is, however, 
not the only difference between pointlike and extended 
T3-scatterers. The FT-patterns also change for small mo- 
menta; most noticeably, the q7-pcaks are enhanced due 
to small-angle scattering of quasiparticles by the smooth 
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FIG. 2: Effect of out-of-plane scattering potentials of different type on the power spectrum. Shown is the power spectrum for: 
1st column: pointlike T3-scatterer with strength V=120meV, 2nd column: extended T3-scatterer with Vb=40meV, A = r z = 1.5a, 
3rd column: pointlike n-scatterer with SA = Ao on the four central bonds, 4th column: slightly extended ri-scatterer with 
8A — Ao, A = r 2 = 1.2a, 5th column: combination of 7.5% extended ri-scatterers (SA = Ao, A = r z = 1.2a) and 3% extended 
T3-scatterers (Vb=40meV, A = r z = 1.5a), 6th column: same as 5th column plus 0.2% unitary (V=3.9meV) pointlike scatterers, 
7th column: experimental power spectrum [l2lll3|l . 
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potential. Whereas the q7-peaks appear only as tiny 
spots in the case of a pointlike impurity pr| . too small 
to be resolved in Fig. they become broad spots in the 
case of an extended impurity and are clearly visible at the 
energies |cj| = 14meV, 22meV in Fig. [2] This supports 
the proposal by Zhu et al. j2|| that the enhancement of 
forward scattering due to spatially extended impurities 
broadens the octet peaks. 

In addition to the existence of weak extended T3- 
scatterers, which could arise e.g. due to random sub- 
stitutions of Sr by Bi, the pair interaction is most likely 
modulated due to the presence of non-stoichiometric oxy- 
gen dopant atoms |3l| . This results in inhomogeneities 
in the magnitude of the superconducting order parame- 
ter and enhanced Andreev scattering which can be mod- 
eled as r\ -scattering within a single- impurity approach. 
The 3rd column in Fig. [5] shows the power spectrum for 
a "pointlike" ri-scatterer, i.e., a modulation of the or- 
der parameter on the four bonds surrounding the im- 
purity site as SA(0,±x) = 8A(±x,0) = -<5A(0,±y) = 
— 5A(±y, 0) = SA. Obviously, the resulting FT-pattern 
looks very different from a pointlike T3-scattercr. The 
background features, i.e., the finite FTDOS at wavevec- 
tors that do not correspond to one of the octet peaks, are 
much less pronounced. At high energies (\u>\ = 30mcV, 
38meV), the power spectrum is clearly dominated by qi- 
peaks, whereas the q7-peaks are completely absent even 
at the lowest energies. This behavior can be understood 
in terms of the momentum dependence of the scatter- 
ing matrix clement for a "pointlikc" ri-scattercr. Fourier 
transformation of the "pointlike" order parameter mod- 
ulation yields 



<5A kk < = <SA(Ak + Ak')/Ao , 



(9) 



i.e., ^Akk' vanishes for wavevectors k and k' which con- 
nect points of the Brillouin zone where A k and A k ' have 
opposite signs but equal magnitudes. Among the octet 
vectors, only qi, q4 and qs are allowed scattering pro- 
cesses contributing to the FTDOS for a "pointlike" t\- 
scatterer, whereas q2, q3, q6 and <\i arc suppressed be- 
cause the matrix element <5Ak.k' vanishes for these q, see 
Fig. ^ The 4th column of Fig. [21 shows the power spec- 
trum for the case of a slightly more extended n-scatterer. 
The primary effect of extending the spatial range of the 
order parameter modulation is analogous to the T3-case, 
i.e., the intensity at large momenta and background fea- 
tures are suppressed resulting in even more pronounced 
qi-peaks. 

The 5th column of Fig. [21 presents the power spectrum 
of an extended impurity with a combined t\- plus 13- 
potcntial, where each single component corresponds to 
the potential used for the extended T3/ri-case in the 2nd 
/4th column. Since we assume that the extended T3- 
scatterers result from the 3% Sr/Bi-disorder and the ex- 
tended Ti-scatterers from the 7.5% dopant oxygens the 
FTDOS of each component has been weighted with a 
prcfactor corresponding to the ratio of these two concen- 
trations. Alternatively, a scatterer with combined T\- and 




FIG. 3: Cut along the diagonal of the Brillouin zone (k x — 
through the power spectrum of a strong pointlike T3-scatterer 
with V=3. 9eV (solid line), a weak extended T3-scatterer with 
Vo=40meV, A = r z =1.5a (dashed line) and a "pointlike" ri- 
scatterer with 5 A = Ao (dotted line) for an energy of (a) 
o>=14meV and (b) u>=30meV. 



T3-character could also arise due to the fact that each oxy- 
gen dopant, besides modulating the pair potential, most 
likely also possesses a small but finite conventional po- 
tential component. The extended t± + r^-case obviously 
combines the main characteristics of the extended n and 
the extended T3-case. At \u>\ — 14meV the power spec- 
trum is dominated by q7-peaks, at \uj\ = 22meV both 
the q7- and the qi-peaks are present, whereas the high- 
est frequencies \ui\ = 30meV and \oj\ = 38meV are clearly 
dominated by q^peaks. 

In addition to out-of-plane impurities, a realistic de- 
scription of the disorder in BSCCO has to take into ac- 
count the in-plane native defects with concentration es- 
timated as 0.2% from STM. These defects are usually 
modeled as pointlike unitary impurities located within 
the Cu02-plancs. The 6th column of Fig. [21 displays the 
power spectrum resulting from the combined effect of the 
extended n- and 73-scatterers, as displayed in the 5th 
column of Fig. [21 plus 0.2% pointlike unitary impurities 
(here V=3.9eV). The existence of native defects within 
the Cu02-planes of BSCCO is only of minor importance 
for the high-frequency power spectrum, because their 
concentration is much smaller than that of the out-of- 
plane defects. At smaller frequencies, however, the strong 
pointlike impurities produce additional large-momentum 
features, especially close to (it, it). This is precisely what 
is observed in experiment: large- momentum octet vec- 
tors are only observed at low frequencies. The reason why 
the large-momentum features appear only at low energies 
can be understood from Fig. [31 where the power spectrum 
along the Brillouin zone diagonal is compared for a single 
unitary pointlikc scatterer with the power spectrum re- 
sulting from a single weak extended scatterer. Note, that 
the intensity of the large-momentum peak, which results 
from the pointlike unitary scatterers, is almost frequency 
independent (see Fig. [3J). However, its weight relative to 
the small-momentum features of FTSTS signal decreases 
for large bias values. This occurs mainly because the 
intensity of the low momentum features, contributed by 
the weak scatterers, grows rapidly with frequency. 
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Overall, the power spectrum calculated within our re- 
alistic disorder model (6th column in Fig. |2J reproduces 
the experimental one (7th column in Fig. EJ very well. 
This allows us to interpret some of the key features of 
the experimental power spectrum in the following way: 

(i) the dominant qi -peaks arise mainly from scattering by 
the order parameter modulations, i.e., by ri-scattcring, 

(ii) the broadness of the q7-peaks, which are quite pro- 
nounced at small energies, is caused by poorly screened 
out-of-plane potential impurities, which act as extended 
T3-scatterers, (iii) the large momentum features at small 
frequencies originate from the in-plane native defects, 
which are modeled as pointlike unitary scatterers. 

IV. CONCLUSIONS 

We have suggested that a realistic model for the out- 
of-plane disorder in BSCCO-2212 should include: (i) 
order parameter inhomogeneities, which we have pre- 
viously attributed to the fact that the oxygen dopant 
atoms modulate the pair interaction locally [3lj and 
which result in enhanced Andreev or Ti-scattcring, (ii) 
extended weak potential or T3-scattercrs, which could 
originate e.g. from the random substitution of Sr by 
Bi and which are spatially extended due to poor screen- 
ing within the Cu02-planes. In combination with the 
0.2% in-plane native defects observed by STM, which 
we model as pointlike unitary potential scatterers, our 
model for the out-of-plane disorder is able to reproduce 
the essential characteristics of the experimental Fourier 
transform patterns in BSCCO. The fact that the out-of- 
plane scatterers are spatially extended explains the ab- 
sence of large-momentum octet peaks and the suppres- 
sion of "background" features at higher energies. The 
large-momentum features at small energies are caused 



by the 0.2% pointlikc strong in-plane scatterers. The q7- 
peaks are broadened considerably due to enhanced for- 
ward scattering caused by the spatially extended poten- 
tial impurities. The t\ -scattering, which arises from the 
order parameter inhomogeneities observed in BSCCO, 
strongly enhances the dispersive qi-peaks especially at 
higher energies. The reason why n -scattering mostly 
contributes to qi -peaks is a simple consequence of the 
c?-wave symmetry of the corresponding scattering matrix 
elements. 

The agreement of the qualitative features of the power 
spectrum predicted by the model introduced in Ref. 
for the out-of-plane O dopants with experiment is fur- 
ther evidence in support of the contention put forward 
there that the O defects in the BSCCO-2212 system are 
actually modulating the pair interaction g locally. Of 
course comparisons of this type cannot directly determine 
the microscopic origin of the modulation of g, but since 
the experimental data with which we compare are from 
slightly overdoped samples, it seems unlikely that signifi- 
cant inhomogeneity is being driven by competition with a 
second order parameter. This therefore suggests that an 
atomic-scale modulation of electronic hopping parame- 
ters or electron-phonon coupling constants is indeed tak- 
ing place in these materials near the dopant atoms. |3lj . 
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